Quantum phase transition to unconventional multi-orbital superfluidity in optical lattices 



Parvis Soltan-PanahiQDirk-Soren LiihrnannQ Julian Struck, Patrick Windpas singer, and Klaus Sengstock 

Institut fur Laser-Physik, Universitdt Hamburg, Luruper Chaussee 149, 22761 Hamburg, Germany 

Orbital physics plays a significant role for a vast number of important phenomena in complex condensed 
matter systems such as high-Tc superconductivity and unconventional magnetism. In contrast, phenomena in 
superfluids - especially in ultracold quantum gases - are commonly well described by the lowest orbital and 
a real order parametei^ . Here, we report on the observation of a novel multi-orbital superfluid phase with a 
complex order parameter in binary spin mixtures. In this unconventional superfluid, the local phase angle of the 
complex order parameter is continuously twisted between neighboring lattice sites. The nature of this twisted 
superfluid quantum phase is an interaction-induced admixture of the p-orbital favored by the graphene-like band 
structure of the hexagonal optical lattice used in the experiment. We observe a second-order quantum phase 
transition between the normal superfluid (NSF) and the twisted superfluid phase (TSF) which is accompanied 
by a symmetry breaking in momentum space. The experimental results are consistent with calculated phase 
diagrams and reveal fundamentally new aspects of orbital superfluidity in quantum gas mixtures. Our studies 
might bridge the gap between conventional superfluidity and complex phenomena of orbital physics. 



The topological properties of graphene and its remarkable 
band structure have recently opened a new field in physics^. 
The linear dispersion relation at the Dirac points proves to 
be a fascinating key aspect of this material as it gives rise to 
phenomena such as quasi-relativistic particles^ and an anoma- 
lous quantum Hall effect^ ^. The possibility to realize hexag- 
onal optical lattices^ permits the emulation of graphene-like 
physics with ultracold atoms^ ^. In particular, loading bosonic 
particles in such lattices renders completely new possibilities 
such as studying next-nearest neighbor processes and tunnel- 
ing blockades in multi-component systems^. In general, op- 
tical lattices have proven to be a versatile tool to simulate 
Hubbard-like systems and actively drive and monitor quantum 
phase transitions^^ The important role of higher orbitals 
has recently been demonstrated for quantum gas mixtures in 



the strongly correlated regime^SEIl. However, for weakly in- 
teracting systems such as superfluids, higher orbitals are gen- 
erally expected to have only marginal effects. So far, orbital 
superfluidity has been observed only in excited states with 
limited lifetimes ^^ ^^ 

Here, we demonstrate the realization of a bosonic super- 
fluid ground-state where higher orbital physics plays a crucial 
role. In conventional superfluids (Fig. la), the local phase an- 
gle <l>(r) of the order parameter 

^{v) = y^^e^^^'^) (1) 

is constant (represented by arrows) and (p{t) can therefore be 
chosen as real. In contrast, the observed TSF ground-state re- 
veals a non-trivial complex superfluid order parameter, where 
the phase factor e**^'^'^^ is continuously twisted in the complex 




FIG. 1: Normal and twisted superfluid phase. The contour plots depict the superfluid order parameter (^(r) = ^n(r)e*^*^'"'* (red lines) 
in the hexagonal lattice. The arrows and color coding represent the local phase angles $(r). a, Normal superfluid: Single-component 
spin systems exhibit a real superfluid order parameter with constant local phases, b, Twisted superfluid: In binary spin mixtures, the 
local phases $(r) vary continuously between the sites of the hexagonal lattice leading to a complex superfluid order parameter for 
both components. This twisted phase is caused by a complex admixture of the p-band. c, The spin-dependent lattice causes a density 
modulation in dependence on the spin state |F, m^). The s- and p-orbital wave functions as well as the spin-dependent potentials are 
shown as a one-dimensional cut along two adjacent lattice sites. 
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FIG. 2: Symmetry breaking in momentum space, a, The normal superfluid phase shows a six-fold-rotational symmetry due to the 
lattice geometry. Averaged TOF images are shown for different lattice depths (indicated in the images), b, For certain binary mixtures, an 
interaction-induced reduced three-fold rotational symmetry is clearly observed for both constituents being separated by a Stern-Gerlach field. 
The triangular pattern is opposite for the two constituents. The six-fold rotational symmetry is restored for large lattice depths, c, Calculated 
momentum distributions for s- and p-orbital wave functions for |F, ttif = 0), where the p-orbital wave function in momentum space has an 
alternating sign in its first-order peaks. A complex superposition of both orbitals leads to the observed TOF pattern for the TSF phase even for 
small admixtures of the p-orbital (shown for p-band admixtures oin^ — \ — Us — 0.05). 



plane (Fig. lb). As we will demonstrate, this unconventional 
behavior arises from a strong interaction-induced coupling of 
s- and p-orbitals. Most strikingly, even a quantum phase tran- 
sition between the NSF and the TSF phase is directly observed 
in our experiment. It is driven by the competition between 
intra- and interspecies interactions in binary mixtures. For its 
realization, we use a mixture of ultracold atoms in two spin 
states |F, mp) of the hyperfine manifold F = 1, 2 with Zee- 
man states vfiF (in our experiment ^^Rb atoms^^, see meth- 
ods). The repulsively interacting atoms are confined in a spin- 
dependent hexagonal optical lattice^. 

As a central aspect, the formation of the twisted superfluid 
phase originates from both the spin-dependency and the spe- 
cific topology of the hexagonal lattice. The topology leads to 
a graphene-like band-structure with the particular feature that 
s- and p-bands are separated only on the order of the tunneling 
energy. In addition, the spin-dependency induces an individ- 
ual sublattice structure for different |F, m^) states. This leads 
to an alternating density modulation for |F, m^? ^ 0) spin 
states (see Fig. Ic), which substantially alters the interspecies 
interactions. The combination of both effects causes a strong 
coupling of s- and p-bands for the case of spin mixtures. 

In the following, we first explain how the TSF phase can be 
identified experimentally and subsequently discuss theoretical 
phase diagrams as well as experimental results in detail. The 
NSF phase possesses the expected six-fold rotational symme- 
try in momentum space. This is observed in experiments via 
time-of-flight (TOF) imaging as exemplarily shown in Fig. 2a 
for the single-component spin- state |1, — 1). In stark contrast, 
the twisted superfluid phase is accompanied by a symmetry 
breaking in momentum space which appears as an alternating 
pattern in the first order momentum peaks (Fig. 2b). This re- 



duced three-fold rotational symmetry reflects the occurrence 
of a twisted complex phase factor e**^^'^^ which we observe 
only for mixtures of two spin states. Figure 2b shows the re- 
sults for a 1 : 1 mixture of 1 1, — 1) and 1 1, +1) atoms, where the 
spin states are separated in the experiment by a Stern-Gerlach 
field. The occurrence of the twisted superfluid phase is clearly 
visible for very low lattice depths. Here, both components 
show a complementary momentum distribution reflecting the 
opposite phase twist as indicated in Fig. lb. For increasing 
lattice depths, the transition to the normal superfluid phase is 
observed by the restoration of the six-fold rotational symme- 
try (Fig. 2b). Finally, the overall interference contrast vanishes 
as the system approaches the Mott insulator phase, where the 
atoms localize on individual lattice sites. 

Remarkably, the clear signature of the TSF phase persists 
even for p-band admixtures as small as a few percent, which 
is a typical value in our experiment. This relies on the fact 
that the first order momentum peaks of the p-orbital are much 
stronger than those of the s-orbital leading to a strong am- 
plification of the p-band contributions, which can be clearly 
identified in the experimental signal (see Fig. 2c). The ob- 
served alternating pattern for the twisted superfluid is caused 
by the quantum interference of s- and p-band contributions, 
which permits an extremely sensitive probing of the local 
phase properties via time-of-flight. 

The connection between the hybridization of s- and p- 
orbitals and the transition to the TSF phase can be explained 
as follows. In general, a superposition between s-orbital \s) 
and p-orbital \p) can be written as 

|(^) = V^|5)+e^^V^|p), (2) 

where the coefficients Us and denote the fraction of atoms 
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FIG. 3: Quantum phase transition to a twisted superfluid phase, a, Zero-temperature phase diagram for the |1, 0) state in admixture 
with the |1, +1) state. The diagram is parameterized by the relative occupation of the p-orbital rip = 1 — ris and the hexagonal lattice depth 
Vq. The white area defines the normal superfluid phase (NSF) and the blue area the twisted superfluid phase (TSF), where the color encoding 
represents the value of the order parameter 0. b, Phase diagram of |1, — 1) in admixture with (A) |1, +1), (B) |2, —2) and (C) |2, +2). For a 
lattice depth of Vo = 5-6£^rec, the system undergoes the transition to a Mott insulator (MI), where the entrance points depend on the particular 
spin mixture^. 



in the s- and the p-orbital, respectively. The global phase an- 
gle 6 between both orbitals is crucial for the formation of the 
twisted superfluid. It takes the value which minimizes the en- 
ergy of the system and can lead to two different physical situa- 
tions: For ^ = (or tt), the system is in the normal superfluid 
phase, where no interference takes place and the alternating 
momentum pattern vanishes. In contrast, 0^0 causes a de- 
structive interference of the first order peaks (Fig. 2c) thereby 
revealing the twist of the local phases $(r). Thus, the global 
phase angle takes the important role of an order parameter 
describing the phase transition between NSF and TSF, where 
the TSF phase is defined by a non-zero value of 0. 

We explore this phase transition theoretically using a multi- 
band mean-field approach (see methods), which leads to the 
phase diagrams presented in Fig. 3. The phase diagrams show 
the results for different binary spin mixtures. In agreement 
with the experimental results, the twisted superfluid emerges 
only in shallow lattices. It is important to mention that in our 
hexagonal optical lattice configuration, different spin states 
preferably occupy different sublattices and therefore the in- 
terplay of intra- and interspecies interaction strongly depends 
on the spin mixture considered. In particular. Fig. 3b demon- 
strates that the TSF phase area is drastically reduced for spin 
mixtures predominately occupying the same sublattice (mix- 
ture C) in comparison to spin mixtures, where each compo- 
nent occupies a different sublattice (mixture A and B). In ad- 
dition, the occurrence of the TSF phase depends on the admix- 
ture of the p-band orbital. This also explains the absence of 
the TSF phase for single-component samples (Fig. 2a) where 
the population of the p-orbital is negligible. 

In the following, we experimentally investigate the phase 
diagrams above. As an experimental indicator characteriz- 
ing the NSF and TSF phase, we define a triangular interfer- 
ence contrast /a which is illustrated in Fig. 4. It serves as 
a measure of the order parameter of the NSF-TSF transi- 
tion, where /a 7^ corresponds to the TSF phase. For spin 
mixture A, the TSF phase is clearly resolved for Vq < 4£^rec 



(Fig. 4a). As expected for symmetry reasons, both compo- 
nents exhibit the same triangular contrast |/a |. In accordance 
with the phase diagrams presented in Fig. 3, mixture B ex- 
hibits a similar behavior as mixture A, where in both cases 
different sublattices are occupied by the constituents (Fig. 4b). 
The substantial difference of the predicted TSF phase areas for 
mixtures B and C is also clearly revealed in our experiment. 

To gain further insight into the underlying processes of the 
NSF-TSF transition, we turn back to its theoretical descrip- 
tion. The observed quantum phase transition is entirely driven 
by the competition between intra- and interspecies interac- 
tions and can thus occur at zero temperatures. We apply mean- 
field theory, where we restrict our analysis to an effective two- 
mode Hamiltonian^^ (see methods). For the superfluid order 
parameter, we consider s- and p-band contributions described 
by Eq. (2), where the order parameter of the transition 6 is the 
relative phase angle between the two orbitals. In the calcula- 
tion, higher bands and non-zero quasi-momentum states can 
be neglected to first order. 

For simplicity, we discuss here a mixture of |1,0) with 
1 1,-1) atoms (see Fig. 3a). In this case, the symmetry of 
s- and p-band wave functions leads to only two competing 
^-dependent terms in the energy functional, namely, the in- 
terspecies interaction Vsp{0) = 2N'^Vsp y/nsTip cos{6) and 
the intraspecies interaction Wsp{0) = 2N'^Wsp UsTip cos(26>). 
Here, N is the number of particles in each component, Vsp 
and Wsp are integrals of s- and p-orbital wave functions with 
Vsp < and Wsp > (see methods). For sufficiently small 
admixtures of the p-band n^, the interspecies interaction Vsp 
dominates and therefore = minimizes the energy, which 
corresponds to a normal superfluid. However, for a critical 
value of Up the energy functional Vsp{0) + Wsp{0) no longer 
exhibits a minimum at 6> = 0. As a central result, this defines 
the phase boundary n^^^ of the NSF-TSF transition, which is 
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FIG. 4: Observation of the NSF-TSF quantum phase transition. 

Absolute value of the triangular interference contrast /a (see inset) 
as a function of the lattice depth Vb for a, spin mixture A and b, spin 
mixtures B and C. The TSF phase is identified by a non-zero value of 
I A , whereas /a vanishes in the NSF phase. The error bars indicate 
the standard deviation of the mean. 



given for the considered case by 




(3) 



The phase transition is of second order as the second deriva- 
tive of the energy functional is discontinuous at this bound- 
ary. In the applied theory, the Hamiltonian is invariant under 
the transformation 0^—0 causing a two-fold degenerate 
ground-state in the twisted phase. However, in the experiment 
always the same of the two ground-states is observed. Fur- 
ther studies are necessary to investigate this non-spontaneous 
symmetry breaking. 

Our study of a new type of complex superfluid phase paves 
the way for further investigations of the interplay between or- 
bital physics and strong correlations. In particular, a possible 
competition between the twisted superfluid and the strongly 
correlated Mott insulator phase can be realized by increas- 
ing the interactions, e.g. by means of Feshbach resonances. 
Moreover, further insight into the roles of intra- and inter- 
species interactions can be gained using binary mixtures con- 
sisting of two different atomic states, where both interactions 
differ considerably from each other. In addition, dynamically 
driven phase transitions may be observable in our systems, by 
preparing a dynamical superposition of s- and p-orbitals for 
one of the two spin components by microwave coupling. 



Methods 

Creation of spin-dependent hexagonal lattices - The spin- 
dependent hexagonal lattice is realized by intersection of three copla- 
nar laser beams under an angle of 120 degrees. The laser beams 
are derived from a Ti: Sapphire laser operated at 830nm (red de- 
tuned), where each beam is linearly polarized within the plane of 
intersection. Orthogonal to the plane, we apply a retro-reflected one- 
dimensional lattice at V\t> — 8.8£^rec operated at the same wave- 
length (for details see Ref. ^ For the experimental parameters, we 
expect an admixture of the p-orbital between zero (for vanishing in- 
teractions) and a few percent for large lattice depths. 

Preparation and detection schemes for spin mixtures - We start 
with a Bose-Einstein condensate of typically several 10^ atoms in 
the stretched state |1, — 1), which is confined in a nearly isotropic 
crossed dipole trap {uj ^ 2ti x 90Hz). The preparation of the dif- 
ferent pure and mixed spin states is performed with aid of radio- 
frequency and/or microwave sweeps. After the state preparation we 
ramp up the optical lattice within 80ms using an exponential ramp. 
Within the ramping time the coherence between different spin states 
is lost. To separate different spin-components during 27ms time-of- 
flight, a Stern-Gerlach gradient field is applied prior to absorption 
imaging. 

Theoretical model - To first order, one spin species experiences 
the interaction with the non-interacting density M\(^' [y)\^ of the 
other spin species with M atoms. Thus, we can write the effective 
Hamiltonian for one spin state as 



+ gnM\^\v)\^ + ^-^^\v)i;{v)U{v) 



(4) 



where Hq — ^ + is the operator for kinetic and potential 
energy and ip is the bosonic field operator. The intra- and interspecies 
interaction parameters are labeled as Qa — ^^^a^ with a = tt? ti^ 
respectively, and ^ lOOao- For shallow lattices, we assume that 
only s- and p-band Bloch functions (ps,p with quasi-momentum q = 
contribute. For large particle numbers and weak interactions, we 
apply mean-field theory and expand the field operators according to 
Eq. (2) 

'0(r) y^(ps + e^^/n^(pp, (5) 

where (ps,p are real functions. The energy functional can be divided 
in a ^-independent and a ^-dependent part, where the latter is given 



by He {rip, 0) 



, + Wsp + Xsp with 



2MNVsp ^/n^cos{0), 
sp risUp cos(2^), 



(6) 



These terms 



4iV (XsUs + XpUp) yyusUp cos(^). 
depend 



Vsp 

Wsr. = 



on 

^p, Xs 



integral 
integrals 



the interspecies 
the intraspecies 

= ^fd^r(pl\(ps\'^^p and 
due 

to reasons of parity. The phase boundary of the phase transition 
between normal and twisted superfluid phase is defined by 



X 



^ fd^r (f*\(pp\^(pp. The latter two vanish for ttif 



2Wsp\/np{l - Up) 



M 
2N 



Vsp + Xs{l- 



+ XpUp 



(J) 



When approaching the Mott insulator transition, higher quasi- 
momentum states become occupied and the two-mode description 
presented here is no longer fully valid. This could explain quantita- 
tive deviations between theory and experiment. The momentum dis- 
tribution yOTOF can be calculated using the Fourier transformed wave 
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